Elliptic Flow from Nonequilibrium Color Glass Condensate Initial Conditions 
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A current goal of relativistic heavy ion collisions experiments is the search for a Color Glass 
Condensate as the limiting state of QCD matter at very high density. In viscous hydrodynamics 
simulations, a standard Glauber initial condition leads to estimate Aurj/s ~ 1, while a Color Glass 
Condensate modeling leads to at least a factor of 2 larger rj/s. Within a kinetic theory approach 
based on a relativistic Boltzmann-like transport simulation, we point out that the out-of-equilibrium 
initial distribution proper of a Color Glass Condensate reduces the efficiency in building-up the 
elliptic flow. Our main result at RHIC energy is that the available data on V2 are in agreement with 
a Airrj/s ~ 1 also for Color Glass Condensate initial conditions, opening the possibility to describe 
self-consistently also higher order flow, otherwise significantly underestimated, and to pursue further 
the search for signatures of the Color Glass Condensate. 
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Ultra-relativistic heavy-ion collisions (uRHICs) at the 
Relativistic Heavy-Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) create a hot and dense system 
of strongly interacting matter. In the last decade it has 
been reached a general consensus that such a state of 
matter is not of hadronic nature and there are several 
signatures that it is a strongly interacting quark-gluon 
plasma (QGP) [H-i]. A main discovery has been that 
the QGP has a very small shear viscosity to density en- 
tropy, "q/s, which is more than one order of magnitude 
smaller than the one of water 0, and close to the 
lower bound of l/47r conjectured for systems at infinite 
strong coupling [6| . A key observable to reach such a con- 
clusion is the so-called elliptic flow, V2 = {cos{2ipp)) = 
{{pI - p1) I ipi + vl)) ■, with ipp being the azimuthal angle 
in the transverse plane and the average meant over the 
particle distribution. In fact, the expansion of the created 
matter generates a large anisotropy of the emitted parti- 
cles that can be primarily measured by V2. Its origin is 
the initial spatial eccentricity, ex — {x^ — y^)/ {x^ -I- y^) , of 
the overlap region in non-central collisions. The observed 
large V2 is considered a signal of a very small 77/ s because 
it means that the system is very efficient in converting ex 
into an anisotropy in the momentum space U2, a mecha- 
nism that would be strongly damped in a system highly 
viscous that dissipates and smooths anisotropies L fl-la l - 
Quantitatively both viscous hydrodynamics [UET 



and transport Boltzmann-like approaches agree 
in indicating an average rj/s of the QGP lying in the 
range Atttj/s ~ 1 — 3. 

Along with the existence of a deconfined QGP mat- 
ter and the understanding of its properties, the uRHIC 
program offers the tantalizing opportunity to explore the 
existence of an exotic state, namely the Color Glass Con- 
densate (CGC) see 0,[2l| for reviews. Such a 
state of matter would be primarily generated by the very 
high density of the gluon parton distribution function at 



low X (parton momentum fraction) , which triggers a sat- 
uration of the gluon distribution function at a pt below 
the saturation scale, Qs 1131 • Even if at first sight sur- 
prisingly, the study of the shear viscosity ij/s oi the QGP 
and the search for the CGC are related. In fact, the main 
source of uncertainty for rj/s comes from the unknown 



initial conditions of the created matter [10|, |23|, |24 1 



A sirnple geometrical description through the Glauber 
model [25| predicts a ex smaller at least 25-30% than the 
eccentricity of the CGC, for most of the centralities of the 
collisions, see for example results within the Kharzeev- 
Levin-Nardi (KLN) model [13, SIH, factorized KLN 
(fKLN) model [28J, Monte Carlo KLN (MC-KLN) model 



and dipole model [2J, [30|. In fact, the satu- 



ration of gluon distribution function is stronger in the 
central overlapping region, making the space distribution 
sharper than the one coming from a simple geometrical 
overlap. Other CGC models like the ones based on the so- 
lution of the classic Yang-Mills equations predict a some- 
what smaller initial eccentricity 31]. The uncertainty 
in the initial condition translates into an uncertainty on 
rj/s of at a least a factor of two as estimated by mean of 
several viscous hydrodynamical approaches [ifl, 0, ■ 
More explicitely, the experimental data of V2{pt) at the 
highest RHIC energy are in agreement with a fluid at 
Airrj/s ~ 1 according to viscous hydrodynamics simula- 
tion, assuming a standard Glauber initial condition. As- 
suming an initial fKLN or MC-KLN space distribution 
the comparison favors a fluid at Anrj/s ^ 2. The reason is 
the larger initial of the fKLN, which leads to larger V2 
unless a large rj/s is considered. However, in [2^ it has 
been shown that viscous hydrodynamics fails to repro- 
duce both V2 and V3 = {cos{3ipp)) if the same Amj/s ~ 2 
is assumed. At variance a Glauber initial condition can 
account for both with the same Atttj/s = 1 

In this Letter, we point out that the implementation 
of the melted CGC in hydrodynamics takes into account 



2 



only the different space distribution respect to a geomet- 
ric Glauber model, discarding the key and more peculiar 
feature of the damping of the distribution for pt below 
the Qs saturation scale. We have found by mean of ki- 
netic theory that this has a pivotal role on the build-up 

of V2- 

We adopt the model which was firstly introduced by 
Kharzeev, Levin and Nardi [26j (KLN model). In partic- 
ular, to prepare the initial conditions of our simulations 
we refer to the factorized-KLN (fKLN) approach as in- 
troduced in H, 2^ , in which the coordinate space distri- 
bution function of gluons arising from the melted CGC 
is assumed to be 



dycPx^ 



J (fpTPA{x±)pB{x±)<i>{pT,x±,y) , (1) 



where $ corresponds to the momentum space distribu- 
tion in the kr factorization hypothesis 32, |33| , 



Att^Nc 1 

^{pT,x±,y) = j^2_\ ^ J (f'krasiQ'^) 

x0b(x2, (pt - feT)^;a;±) . (2) 

Here xi.2 = exp(±j/)/^/s and the ultraviolet cutoff 
Pt = 3 GeV/c assumed in the pr integral in Eq. ((IJ; 
as denotes the strong coupling constant, which is com- 
puted at the scale — max(fe|,, {px — kxY) according 
to the one-loop j3 function but frozen at = 0.5 in the 
infrared region as in [i^,!!^,!!^- In Eq. ^ pa,b denote 
the probability to find one nucleon at a given transverse 
coordinate, pa{x±_) = 1 — {1 — ainTA{x±) / A)^ , where 
Gin is the inelastic cross section and Ta corresponds to 
the usual thickness function of the Glauber model. 

The main ingredient to specify in Eq. ^ is the unin- 
tegrated gluon distribution function (uGDF) for partons 
coming from nucleus A, which is assumed to be: 



i(a;i,fc|;a;j 



0{Qs 



9{kT - Qs) 



A2 



(3) 



where we see the peculiar feature of the CGC that is 
the saturation of the distribution for pr < Q s', sl similar 
equation holds for partons belonging to nucleus B. Fol- 
lowing [2^ we take the saturation scale for the nucleus 
A as 



[X,Xl 



2GeV^ 



Ta{x±) 
1.53pa{x±) 



0.01 



(4) 



with A — 0.28, and similarly nucleus B. This choice 
is the one adopted in fKLN or MC-KLN and in hydro 
simulations [ifl, [111 to study the dependence of V2ipT) 
on rj/s. Using Eqs. (g]) and (P) we find that (Qs) « 1.45 
GcV in the case & = 0, at y = 0, = 200yl GeV 
and X = 0.01, where the average is understood in the 
transverse plane. 
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FIG. 1: Momentum space distributions in a cylinder of radius 
r = 3fm and space rapidity |?7s| < 1 for the case of fKLN 
model at the initial time (thick solid line) compared to the 
the thermal one (thin solid and dashed lines) at their initial 
evolution time tq and to the one of fKLN after a dynamical 
evolution Ar = 0.8fm/c. Inset: time evolution of T • r^^'^ 
for the three different initial conditions considered. All the 
panels refer to Au-Au collisions at ^/s — 200 GeV, with an 
impact parameter b = 7.5 fm. 



We employ transport theory as a base of a simula- 
tion code of the fireball expansion created in relativistic 
heavy-ion collision 16, 17, 35, 36|, therefore the time 
evolution of the gluons distribution function f{x,p,t) 
evolves according to the Boltzmann equation: 

P^SVl = j dT2dTydT2'{h'h'~hh) 

yc\M\''5\pi+p2~Pv -P2') , (5) 

where d^pk — 2Ek{2nYdTk and corresponds to the 
transition amplitude. 

At variance with the standard use of transport theory, 
we have developed an approach that, instead of focus- 
ing on specific microscopic calculations or modelings for 
the scattering matrix, fixes the total cross section in or- 
der to have the wanted 77/ s. To this end one needs an 
analytical relation between 77, temperature, cross section 
and density; as shown in (36l . w\ , the Chapmann-Enskog 
approximation supplies such a relation with quite good 
approximation, in agreement with the results obtained 
using the Green Kubo formula. Therefore, we fix 77/s 
and compute the pertinent total cross section by mean 
of the relation 



Ctot 



{P) 



1 



(6) 



pg{a) i]/s ' 

which is valid for a generic differential cross section 



da/dt ~ a1/{t — m^,)^ as proved in [37[. In the above 
equation a = T/mD^ with to^i the screening mass regu- 
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lating the angular dependence of the cross section, while 



.9(a) = Yo ' '^yy 



2/- + - ) K3{2y) - yK2{2y) 



Hay) , 



with Kn the Bessel function and h corresponding to the 
ratio of the transport and the total cross section. The 
maximum value of (7, namely ^(m^) — >■ 00) = 10, is 
reached for isotropic cross section; a smaller value of g{a) 
means that a higher atot is needed to reproduce the same 
value of r]/s. However, we notice that in the regime were 
viscous hydrodynamic applies (not too large rj/s andpr) 
the specific microscopic detail of the cross section is ir- 
relevant and our approach is an effective way to employ 
transport theory to simulate a fluid at a given rj/s. From 
the operative point of view, keeping ry/s constant in our 
simulations is achieved by evaluating locally in space and 
time the strength of the cross section by means of Eq. ^ , 
where both parton densities and temperature are com- 
puted locally in each cell. 

In the following, we will consider three different types 
of initial distribution function in the phase-space, two of 
which are the one employed till now for the investigation 
of the rj/s, while the third one is the genuine novelty of 
the present study. Furthermore, we refer to Au + Au 
collision at ^/s — 2QQ AGeV and b = 7.5 fm, belonging 
to the 20 — 30% centrality class. In this case, our result 
for initial eccentricity in the fKLN model is tx = 0.357 
(which is in agreement with MC-KLN 28| result used 
in hydro simulations). The standard initial condition for 
simulations of the plasma fireball created at RHIC is a 
a?-space distribution given by the Glauber model and a 
p-space thermalized spectrum in the transverse plane at 
a time r = 0.6 fm/c with a maximum initial temperature 
To = 340MeV. In this case, for a standard mixture of 
Npart and Ncoii we find = 0.284. We will refer to this 
case as Th- Glauber. Instead the study of the impact of 
an initial CGC state has been performed considering an 
a;-space distribution given by the fKLN (or MC-KLN), 
while in the momentum space the spectrum has been con- 
sidered thermalized at tq ~ 0.6fm/c; we refer to this case 
as Th-fKLN and it is represents the case implemented in 
hydrodynamics, that has lead to the conclusion that the 
CGC suggests a 47r?7/s - 2 [l3,[2i|2^. The third initial 
conditions is the full fKLN initial conditions where, be- 
yond the X- space, the saturated distribution in p-space is 
implemented as well, see Fig. [Ijsolid thick line. As initial 
time we take tq = O.lfm/c because in this case there is 
no pre-assumption of thermalization. This is not usually 
considered in hydrodynamics because there it is implic- 
itly assumed a distribution function in p-space in local 
equilibrium, at least in the transverse plane. For all the 
previous cases, as usually done, a Bjorken scaling at the 
initial time is assumed, identifying momentum rapidity 
yz with space rapidity rjs = arctg~^ {z / t) . The particle 



multiplicity, dN/dy, is fixed to correctly reproduce the 
experimental one for all the three cases. 

In Fig. [T] we plot the initial spectra for the fKLN 
(thick solid hue), Th-fKLN (dashed hue) and Th-Glauber 
(7) (thin solid line) at their respective initial times tq, and 
the spectrum of the fKLN model after a time evolution 
At = 0.8 fm/c (dashed green line). We notice that ini- 
tially the fKLN spectrum is quite far from a thermalized 
spectrum; in fact, it embeds the saturation effects which 
are proper of the melted CGC. Nevertheless, the spec- 
trum evolves to a thermalized one within 1 fm/c. Such 
a feature is confirmed by the inset of Fig. [1] where the 
quantity T* ■ t^I^ is shown, with T* = E/iN represent- 
ing, the temperature in the case of a thermalized system. 
It is known that in the case of ID expansion a thermal- 
ized system should keep such a ratio constant. We find 
that in the case of the fKLN (solid green line) the product 
7^*.T-i/3 ig strongly dependent on time because the system 
is quite far from equilibrium; however at r ^ 0.8 fm/c 
both the value and the time evolution become indistin- 
guishable from the thermal cases represented by the Th- 
Glauber and Th-fKLN. We notice a little adjustment also 
for these cases that we have indicated as thermal. The 
reason is that the initial spectra are thermal only in the 
transverse plane, but there is a boost invariance in the 
longitudinal direction. This causes a little re-adjustment 
that would disappear assuming a thermal spectral also 
in the longitudinal direction. 

Our results on thermalization time are not in disagree- 
ment with earlier studies showing that two-body colli- 
sions are insufficient to achieve a fast thermal equilibrium 
(isj . In fact, in that case a perturbative QCD two-body 
cross section is employed, corresponding to rj/s about 
one order of magnitude larger than in our case. The dif- 
ference is of course that we normalize the cross section 
to get the wanted 7//s, and our scattering rates even if 
obtained effectively by two-body scattering is very large, 
leading eventually to a fast thermalization. 

In the left panel of Fig. [2] we plot the V2{pt) for the 
case of Th-Glauber (thin solid line) and Th-FKLN (thick 
solid line) at a fixed Anrj/s = 1. The Glauber initial 
condition reproduces quite well the data (circles); in the 
case of Th-fKLN (thick solid line) one gets a too large 
f 2 and for such initial conditions the agreement with the 
data is achieved only if the rj/s is increased by a factor of 
two (dashed line). These results are in agreement with 
the ones obtained from viscous hydrodynamics 13, 23, 
[j^ l , showing the solidity and consistency of our transport 
approach at fixed 77/ s. 

In the right panel of Fig. [2] we present our novel result 
for the fKLN model, when the CGC distribution function 
is implemented in both the x and p spaces. We find that 
fKLN with a Anrj/s — 1 (thick solid line) gives a V2{pt) 
quite similar to the Th-Glauber, while in such a case if 
Airq/s = 2, dashed line, the V2{pt) would be too small. 
Our interpretation is that the initial larger is compen- 
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FIG. 2: Elliptic flow V2{pt) for the different initial conditions 
and rj/s as in the legend. All the calculations refer to Au-Au 
collisions at ^/s — 200 GeV, with an impact parameter b = 7.5 
fm. 




T-Tg [fm/c] 

FIG. 3: Evolution of V2 at pt ~ 1.5 GeV as a function of the 
evolution time for all the different initial conditions and i]/s 
values. Calculations refer to Au-Au collisions at = 200 
GeV, with an impact parameter b = 7.5 fm. 



sated by the key feature of an almost saturated initial 
distribution in p-space below the saturation scale Qs- In 
fact, we have already checked that increasing artificially 
the Qs results in a further lowering of V2{pt) (details will 
appear in a forthcoming publication). In other words the 
initial out-of-equilibrium fKLN distribution reduces the 
efficiency in converting Cx into V2. In fact, the elliptic 
flow can be understood as a larger slope of the momen- 
tum spectrum in the out of plane x direction respect to 
the y one caused by a larger pressure in the x direction 
due the elliptical shape. The net effect in terms of the 
difference of the particle yields between the two direc- 
tions is larger if the spectra are decreasing exponentially 
respect to the case in which they are nearly flat as a 
function of pt- A detailed study is in preparation, but 
the result we present in this Letter shows that the initial 
out-of-equilibrium function implied by the CGC cannot 
be discarded in studying the build of the V2 and the ef- 
fect change significantly the estimate of the rj/s oi the 
system. 

In Fig. [21 we show the build-up of V2 normalized to 
their respective asymptotic value (r oo) at a fixed 
Pt = 1.5 GeV as a function of time for all the initial con- 
ditions considered. The behavior splits into two classes: 
one is represented by the thermal distributions; the other 
one is given by the case in which the initial spectrum is 
the fKLN one. More specifically, we find that the mag- 
nitude of the rate of increase of V2 for the thermal distri- 
butions is large from the very beginning of the evolution, 
while it is quite reduced for fKLN. We further note that 
after about 1 fm/c roughly corresponding to the thermal- 
ization also for the full fKLN distribution the increase of 
V2 with time becomes very similar to Th-fKLN and Th- 
Glauber. This observation further confirms that it is the 
initial out-of -equilibrium and nearly saturated distribu- 
tion that dampens the efficiency in converting the space 
eccentricity into the V2. Therefore, even if the thermaliza- 



tion sets in quite quickly, as assumed in hydrodynamics 
T ^ 0.8 fm/c, such a time duration cannot be neglected 
in studying the impact of initial conditions, and the ex- 
istence of a CGC state. 

In conclusion, our study based on kinetic theory shows 
that the elliptic flow in a system depends not only on 
the pressure gradients and the rj/s of the system, but 
also on the initial distribution in momentum space. A p- 
distribution with a saturation behavior generates smaller 
V2 respect to the thermal one. This result is quite gen- 
eral, and we expect it should be valid, besides QGP in 
uRHICs, for systems like cold atoms in a magnetic trap 
which are characterized by a value of r]/s close to the 
QGP one [s^. In the specific case of the CGC mat- 
ter studied here, the effect of the initial nonequilibrium 
distribution affects the estimate oi rj/s of about a fac- 
tor of two. The relevance of our results is further en- 
hanced by the fact that Th-fKLN with iirri/s ^ 2 would 
generate a low respect to the available data. This 
is the main conclusion of 23|, namely that Th-fKLN is 
not able to describe for the experimental observations, 
might be revised. Our results open the way to study 
the effect of p-distribution functions on higher harmon- 
ics Vn = {cos{nipp)) as a promising tool to further nar- 
row the evaluation of the ij/s of the QGP, similarly to 
the study of anisotropics in the cosmic microwave back- 
ground (CMB). 
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